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Abstract Creations of light anti-nuclei (anti-deuterium, anti-tritium, anti- 
{3JQ[ He3 and anti-He4) are observed by collaborations at the LHC and RHIC ac- 

celerators. Some cosmic ray experiments are aimed to find the anti-nuclei in 
cosmic rays. To support the experimental studies of the anti-nuclei a Monte 
Carlo simulation of anti-nuclei interactions with matter is implemented in the 
Geant4 toolkit. The implementation combines practically all known theoret- 
ical approaches to the problem of antinucleon-nucleon interactions. 
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1 Introduction 

One of the most exciting puzzles in cosmology is connected with the question 
of the existence of anti-matter in the Universe. A number of dedicated cosmic 
ray experiments aim to search for anti- nuclei [U[2,.3,4 . Also, anti- nuclei have 
been observed in nucleus-nucleus and proton-proton collisions by experiments 
at the RHIC [MMZ] and LHC accelerators [8]. The STAR collaboration at 
RHIC reported in March 2011 that the anti-He4 nuclei were identified in high 
energy nucleus- nucleus collisions (see [5]). The ALICE collaboration at LHC 
confirmed the STAR results in May 2011. 
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An experimental study of anti-nuclei requires a knowledge of anti-nucleus 
interaction cross sections with matter. The cross sections are needed to esti- 
mate various experimental corrections, especially those due to particle losses 
which reduce the detected rate. In practice, various phenomenological ap- 
proaches are applied in order to estimate the antinucleus-nucleus cross sections. 
Thus, a first task is a creation of reliable estimations of the cross sections. Here 
we use the Glauber approach. 

It is obvious that an annihilation can take place at an interaction of an 
anti-nucleus with a nucleus. A lot of mesons can be produced in this way. 
Thus, we have to simulate the meson production. We do this in the framework 
of the quark-gluon string model. 

Low energy mesons can have secondary interactions in nuclear residues. 
We take them into account using the binary cascade model of the Geant4 

toolkit pro]. 



2 Antinucleus-nucleus cross sections 

Anti-proton elastic scattering by deuterons was considered in the classic paper 
by V. Franco and R.J. Glauber [TT]. O.D. Dalkarov and VA. Karmanov [12] 
showed that elastic and inelastic (with excitation of nuclear levels) anti-proton 
scattering by C, Ca, and Pb nuclei are described quite well at p kinetic energies 
of 46.8 and 179.7 MeV. The first calculations of the cross sections of anti- 
deuteron interactions with nuclei in the eikonal approximation were presented 
by Buck et al. [T3] (see also [2]). Cross sections of antideuteron-deuteron 
interactions at Pj = 12.2 GeV/c were calculated using the Glauber approach 
in |15j . They were in good agreement with the experimental data. We use the 
Glauber approach to calculate the antinucleus-nucleus cross sections. 

The elastic scattering amplitude of an anti-nucleus containing A anti- 
baryons on a target nucleus with mass number A is given as |16) : 

- ( A A 

[ <=ij=i 

(1) 

bPjuQ) Mqb)db, 





where the same nomenclature as in [161112] is used. 

The main ingredients of the approach are parameterizations of energy 
dependencies of total and elastic antinucleon-nucleon scatterings. We have 
considered the question in our recent publication |17j . Using them, one can 
calculate various cross sections of antiproton-nucleus and antinucleus-nucleus 
interactions. 

We show in Fig.[T]our calculations in a comparison with experimental data. 
For projectile anti-deuterons we present data at two energies ([H]) ~~ 13.3 and 
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Fig. 1 Absorption cross sections of anti-proton and anti-deuteron interactions with nuclei. 
The points are experimental data |18| . the lines are our calculations. 



25 GeV/c (open and close circles, correspondingly) and calculation results 
at the energies (solid and dashed lines). As seen, the agreement between the 
experimental data and the calculations is rather good. 

The profile function of the elastic antinucleus- nucleus scattering, Paa^), 
is very close to the 2-dimensional Fermi- function, Paa{^) = V[l + exp((b — 
R) /c)} . The Fourier-Bessel transform of the function gives the elastic scattering 
amplitude in the momentum representation. We use analytical expressions for 
an expansion of the transform presented in f |19|) to calculate differential elastic 
antinucleus-nucleus cross sections. 

Calculations of the total, inelastic and elastic cross-sections of of an anti- 
nucleus interactions with nuclei allow to be determined a point where the 
anti-nucleus penetrating through a matter will interact. 



3 Simulation of multiparticle production in antinucleon-nucleon 
interactions 

The main channel of the antinucleon-nucleon interactions at low energies is 
the annihilation into 3, 4, or 5 mesons. It is commonly assumed that the 
reaction is going through the re-arangement of quarks and anti-quarks in the 
colliding particles, see Fig. [5^. At high energies, the reaction will result in a 
three quark-gluon string creation. 

If a quark and an antiquark annihilate (Fig.[2jD), diquark-antidiquark string 
will be left. Such a string must fragment into a final state containing a baryon 
and an anti-baryon. A description of the fragmentation of such strings is one 
of the problems in high energy physics. We have solved it using experimental 
data on the reactions pp — > nn, AA, finir , and so on. 

A string junction couples quarks in a baryon. String junction's annihilation 
takes place in the process of Fig. which results in a three string creation. 

At annihilation of quark, anti-quark and string junctions, two strings will 
be created (see Fig. [3b). If the energy is sufficiently low, two mesons will be 
produced. 
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Fig. 2 Quark flow diagram. 



The diagrams " d" and " e" describe the processes with one string creation. 
The ordinary quark- antiquark string can fragment at least into two mesons, 
and it can be responsible for the reactions with two final state mesons. A 
hybrid meson with hidden baryon number can be produced in process " d" at 
low energies. 

It is needed to add the " standard" Fritiof model [20] diagrams " g" at high 
energies. One or two strings can be created in the corresponding processes 
without baryon number exchanges. 

The main issue of the considered approach is the energy dependencies of the 
process cross sections. We did not find an acceptable solution of the question. 
Thus, we parameterized the cross sections as follow: 



16 



vs — 4m 2 



[(«- 4m 2 )-°' 175 + 3.125 * (1 - 1.88/Vi)] (mb) (2) 



a h = 3.13 + 140 -(Mth-yG) 2 - 5 (mb), y/s < M th = 2.172 (GeV) (3) 
a h = 6.8/Vs (mb), y/s > M th = 2.172 (GeV) 

a e = 23.3/s (mb) (4) 

a c = a d = (jf =0 (5) 

a g = 35 * (1. - 2.1/Vs) (mb) (6) 

The parameterizations of the cross sections and LUND string fragmenta- 
tion algorithm are implemented in FTF generator of the Geant4 toolkit [10] . 
Using this generator, we obtain results presented in Fig. [31 As seen, the ex- 
perimental dependencies are reproduced, though a more refined consideration 
is needed. 
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4 Antiproton-nucleus interactions 

Cross sections (cr„) of various processes with different multiplicities (v) of in- 
volved nuclear nucleons can be determined using the asymptotic Abramovsky- 
Gribov-Kancheli (AGK) cutting rules at high energies in application to elastic 
scattering amplitude: 



A 

a pA = a ^ (T » = C A d 2 b -r / g{b - s)p A (s, z)d 2 s dz 

v=\ J I J 

1 — J g(b — s)pa(s, z)d 2 s dz 



(J) 



5 (b) = 7 (b) + 7 *(b)- 7 (b) - 7 *(b). 



At low energies, we have introduced finite energy corrections to the AGK- 
cutting rules and obtained 



°M=E°£' "™» = bWPo] + l, p Q ~2(GeV/c). (8) 



d 2 bl 1 



1 — — I g(b — s)pa{s, z)d 2 s dz 



(9) 



1 — I g(h — s)pa(s, z)d 2 s dz 



-v)A/v„ 



Thus at a low energy, only one inelastic interaction of a projectile in a 
nucleus can happen, but there can be other interactions caused by secondary 
particles. The secondary particle interactions are mainly responsible for slow 
neutron and proton productions (see Fig. 2]). We have checked this using two 
variants of the model - FTF with the binary cascade model (BIC) and FTF 



6 



A. Galoyan, V. Uzhinsky 




with the precompount de-excitation model (PRECO) [33] ■ A secondary parti- 
cle cascading and a nuclear residual de-excitation are considered in the binary 
cascade model, the second model does not take into account the cascading. 
Both models BIC and PRECO are presented in the Geant4 toolkit. As seen 
from a comparison of experimental data |21U22j to the model calculations in 
Fig. |4j the finite energy corrections and the secondary interactions are very 
important for understanding experimental regularities. 



5 Antinucleus-nucleus interactions 

To simulate the distribution (8), we use the following algorithm: starting with 
the expression (7) we ascribe a projectile a power P = v max . A probability 
of an interaction with the first nucleon is equal P jv max . The power decreases 
after the interaction on unit. The probability of an interaction with the second 
nucleon is equal to Pjv maxi where P = v max — 1. If the second interaction 
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has happened, the power is decreased one more. In other case, it is left on the 
same level. This is applied for each possible interaction. After an annihilation 
P is set to zero. One can check that this leads to Exp. (9). 

The same algorithm is applied in the case of antinucleus-nucleus interac- 
tions using the Glauber inelastic cross sections analogous to Eq. (7). A cas- 
cading of secondary particles in the light projectile anti-nuclei is neglected. A 
result of calculations is presented in Fig. [5j where the points are experimental 
data [55], the line is our calculation. 




Fig. 5 Charged particle multiplicity distribution in d + Ta interactions at = 6.1 GeV/c. 



There is no other data on light anti-nucleus interactions with heavy mate- 
rials except [26], though there are a lot of detailed data on dd- interactions, 
and only few data on general properties of the interactions. 



Conclusion 

Using the described approach we have developed in the Geant4 framework 
a Monte Carlo model for simulation of antinucleus-nucleus interactions for 
the projectiles p, d, i, 3 He, He. The model is valid between 100 MeV/c and 1 
TeV/c per antinucleon. A comparison of the model calculations to available 
data shows a good agreement sufficient for most applications in cosmic ray 
experiments and in large HEP experiments including those at the LHC and 
RHIC 
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The Geant4 toolkit is now able to simulate the ant inucleus- nucleus inter- 
actions for all target nuclei since version 9.4-p01. 
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